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摘 要 : 基于 多 种 遥感 数据 ,运用 Mann-Kendall 法 .Theil-Sen median 趋势 分 析 .Pearson 相关 分 析 水 量 平衡 原理 ， 


探讨 2000 一 2020 年 伊犁 河 一 巴尔 喀什 湖 流域 实际 蒸 散发 时 空 变化 规律 及 其 主要 影响 因素 ,并 讨论 了 流域 生态 系统 
水 资源 供给 量 的 变化 。 结 果 表 明 :(1) 流域 上 .中 下游 的 多 年 平均 年 蒸 散 发 量 分 别 为 439.0 mm 317.9 mm 201.1 
mm; HP E 中游 在 夏季 的 日 蒸 散 发 量 最 大 ,而 下 游 在 春季 最 大 ;流域 上 .中 游 燕 散发 量 的 年 内 分 配 均 为 " 单 峰 型 ”， 
峰值 分 别 在 7 月 与 6 月 ,下 游 为 “ 双 峰 型 ”, 峰 值 分 别 在 3 月 与 11 月,(2) 流域 上 .下 游 的 年 蒸 散发 量 均 呈 现 显著 的 上 
升 趋势 ,显著 区 域 主要 分 布 在 上 游 的 伊犁 河谷 、 天 山 山脉 与 下 游 的 伊犁 河 三 角 洲 附近 ;相对 2000 一 2010 年 ,2010 一 
2020 年 伊犁 河谷 和 伊犁 河 三 角 洲 地 区 多 年 平均 年 蒸 散 发 增加 超过 10%。(3) 流域 上 .中 游 蒸 散发 与 气温 和 NDVI 呈 
现 较 高 的 正 相 关 性 ;流域 下 游 蒸 散发 与 土壤 湿度 呈现 较 高 的 正 相 关 性 ,(4) 流域 生态 系统 的 水 资源 供给 服务 总 量 在 
2000 一 2020 年 有 所 下 降 , 其 中 上 游 已 多 次 出 现 缺 口 , 中 、 下 游 在 2020 年 开始 出 现 缺 口 , 需 通过 控制 用 水 总 量 、 提 升 用 


水 效率 保障 水 资源 供需 平衡 。 


关键 词 : 伊犁 河 一 巴尔 喀什 湖 流 域 ; 蒸 散 发 ;时空 变 化 ; 驱动 因素 ; Pearson 相关 


实际 蒸 散 发 (Evapotranspiration ,ET) 包 括 土壤 
藻 散 发 .植物 截留 水 分 蒸发 和 植被 蒸腾 ,是 陆地 水 


流域 的 蒸 散发 呈 下 降 趋势 , 且 同 等 降水 增幅 条 件 
下 ,塔里木 河流 域 的 降水 对 莹 散发 贡献 更 大 ;刘洋 


文 循环 和 能 量 循 环 的 重要 环节 。 位 于 干旱 半 干 旱 
地 区 的 伊犁 河 一 巴尔 喀什 湖 流 域 ( 简 称 “ 伊 巴 流 
域 ”) 是 我 国 与 哈萨克 斯 坦 之 间 的 重要 跨 界 流域 , 流 
域 面积 大 且 流 域内 地 貌 气 候 特征 差异 悬殊 。 随 着 
近年 来 巴尔 喀什 湖 (简称 “ 巴 湖 ”) 水 位 波动 ,伊犁 
河 的 三 角 洲 生态 环境 变化 ,水 资源 分 配 问 题 变 得 
极为 敏感 。 在 全 球 气 候 变 化 和 人 类 活动 加 剧 的 情 
BLP , 蒸 散 发 作为 干旱 半 干 旱地 区 水 量 平 衡 中 最 大 
的 支出 项 ,其 时 空 演变 及 驱动 因素 对 伊 巴 流域 的 水 
资源 分 配 与 管理 意义 重大 。 

已 有 许多 学 者 针对 干旱 半 和 干旱 区 蒸 散发 时 空 
变化 及 驱动 因素 开展 研究 。 李 修 仓 "通过 研究 3 个 
不 同 气候 区 的 典型 流域 发 现 ,1961 一 2010 年 , 塔 里 
木 河 流域 的 实际 蒸 散 发 呈 增 加 趋势 ,而 珠江 和 海河 


收 稿 日 期 : 2022-04-18; 修订 日 期 : 2022-06-25 


等 "利用 相关 性 分 析 得 出 ,中国 西 北 干 旱 区 1960 一 
2019 年 的 实际 蒸 散 发 呈现 增加 趋势 , 且 主 要 受气 
温 、. 相 对 湿度 和 风速 的 影响 ;Xia 等 通过 对 东亚 草 
原生 态 系统 的 研究 得 出 ,内 蒙古 锡林郭勒 地 区 的 落 
散发 和 降水 相关 性 很 高 ;Zhang 等 ”人 研究 发 现 ,内 蒙 
古 中 部 和 西部 地 区 的 蒸 散 发 变化 主要 受 相 对 湿度 
驱动 ; 闫 俊杰 等 "研究 2001 一 2015 年 伊犁 河谷 草地 
蒸 散发 变化 发 现 , 多 年 平均 ET 为 395.74 mm; 植 被 
覆盖 度 、 降 水 及 气温 的 降低 均 为 草地 ET 降低 的 驱动 
因素 ; 邓 兴 泗 等 中 研究 2000 一 2014 年 天 山地 区 燕 散 
发 发 现 ,天 山 山 区 ET 值 的 减 小 趋势 是 由 区 域 降 水 量 
减少 所 致 ; 梁 红 闪 等 3 研究 2000 一 2014 年 伊 巴 流域 
的 蔡 散 量 发 现 ,流域 实际 蒸 散发 分 布 不 均 ,境内 部 
分 多 年 平均 年 ET 为 401.18 mm, 远 高 于 境外 的 
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194.45 mm ,供水 不 足 是 制约 流域 植被 蒸 散 最 主要 的 

干旱 半 干 旱 区 不 同 流域 的 蒸 散 发 变化 受 不 同 
因素 驱动 , 且 前 人 在 伊 巴 流域 的 燕 散 发 研究 方面 已 
取得 相关 研究 成 果 ,但 由 于 数据 源 的 不 同 .时 间 空 
间 斥 度 的 不 同 ,尚未 系统 揭示 伊 巴 流域 的 蒸 散发 演 
变 特征 及 影响 因子 变化 特征 。 因 此 ,本 文 依据 巴 湖 
的 主要 水 量 来 源 一 一 伊犁 河 来 划分 伊 巴 流域 的 上 、 
中 下游, 基于 长 序列 的 遥感 资料 ,揭示 环境 变化 下 
伊 巴 流域 实际 蒸 散发 的 演变 规律 ,并 利用 相关 性 分 
析 谢 析 多 重 环境 因子 对 其 产生 的 影响 ,结果 可 为 伊 
巴 流域 的 水 资源 规划 管理 及 中 哈 跨 界 水 资源 合理 
分 配 提 供 参 考 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

伊 巴 流域 地 理 位 置 见 图 1。 巴 湖 位 于 哈萨克 斯 
坦 的 内 陆 湖 ,伊犁 河 的 径流 是 其 水 量 最 重要 的 补给 
来 源 , 占 总 和 人 湖水 量 的 78.4%51。 伊 犁 河上 游 的 天 
山地 区 是 产 流 区 ,下游 哈萨克 斯 坦 境内 则 是 径流 的 
消散 区 。 以 伊犁 河 三 道 河 子 水 文 站 和 哈萨克 斯 坦 
的 伊犁 村 ( 即 现在 的 卡 普 恰 盖 镇 水 文 站 ) 为 分 界 点 ， 
将 伊犁 河 一 巴尔 喀什 湖 流 域 分 为 上 .中 下游 三 部 
分 ,其 中 下 游 还 包括 除 伊犁 河 外 的 卡拉 塔 尔 河 、 阿 
克 苏 河 . 列 普 瑟 河 、 阿 亚 古 兹 河 的 集 水 面积 中 。 
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流域 高 程 土地 利用 、 植 被 归 一 化 指数 (Normal- 
ized Difference Vegetation Index, NDVI) 降雨 的 空间 
分 布 见 图 2。 由 图 2 可 见 ,流域 高 程 落差 极 大 , 且 降 
雨量 从 山区 向 平原 减少 现象 极为 明显 ,高 山区 多 年 
平均 降水 量 大 于 1000 mm, 低 山区 年 降水 量 为 300~ 
400 mm, 巴 湖 沿岸 年 降水 量 仅 为 150 mm 左右 ”。 
NDVI 的 高 值 区 主要 位 于 流域 上 游 .中 游 高 海拔 地 
PX BEAR URS FRAUD ,而 巴尔 喀什 湖 周边 除 伊犁 河 三 
角 洲 外 植被 覆盖 度 极 低 , 土 地 利用 的 空间 分 布 也 印 
证 了 NDVI 的 分 布 特征 。 
1.2 数据 来 源 
1.2.1 ET、NDVI 数 据 首先 ,采用 美国 航空 航天 局 
(NASA) 提 供 的 APPEEARS 工具 (https:/Npdaacsve. 
cr.usgs.gov/appeears/) 下 载 编号 为 MOD16A2CF (时 
空 分辩 率 为 8 d、500 m) 的 蔡 散 发 数据 和 编号 为 
MOD13C2( 时 空 分 辩 率 为 月 、0.05°) 的 NDVI 数 据 ; 
然后 ,经 过 去 除 无 效 值 重 采样 .人 研究 区 域 提取 和 数 
据 合 并 等 预 处 理 后 ,得 到 栅 格 大 小 为 0.1° 的 月 平均 
数据 ;最 后 通过 Arcmap 软件 将 各 月 的 蒸 散 发 数据 进 
行 逐 像 元 计算 ,分别 得 到 年 与 季节 尺度 的 ET。 
1.2.2 降雨 气温、 相对 湿度 、 风速 数据 由 NASA 
提供 的 饥荒 预警 系统 网 络 土 地 数据 同化 系统 (The 
Famine Early Warning Systems Network Land Data As- 
similation System, FLDAS) 数据 集 , 提 取 研 究 区 
2000 一 2020 年 空间 分 辨 率 为 0.1 的 月 降水 .气温 、 相 
对 湿度 .风速 数据 。 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2021)5452 号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 伊 巴 流域 地 理 位 置 
Fig. 1 Location of the Ili River-Balkhash Lake Basin 
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图 2 伊 巴 流域 基础 信息 空间 分 布 
Fig. 2 Spatial distribution of basic information in the Ili River-Balkhash Lake Basin 


1.2.3 土壤 湿度 数据 NASA 发 布 的 全 球 陆 地 数据 
同化 系统 (The Global Land Data Assimilation Sys- 
tem,GLDAS) 提 供 土 层 厚度 分 别 为 0~10 cm 10~40 
cm、40~100 cm 和 100~200 cm 的 土壤 湿度 数据 ,为 
了 反映 影响 ET 的 供水 条 件 ,同时 考虑 植被 根系 的 生 
长 深度 , 故 从 GLDAS 中 提取 空间 分 辩 紊 为 0.25° 的 
月 尺度 土壤 湿度 (10~40 ecm) 数据 ,经 过 Arcmap 软件 
重 采 样 得 到 空间 分 辨 率 为 0.1° 的 月 土壤 湿度 数据 。 
1.2.4 土地 利用 数据 ”土地 覆 被 数据 选用 欧洲 航天 
局 (European Space Agency, ESA) 气 候 变 化 倡议 
(Climate Change Initiative, CCI) 与 哥 白 尼 和 气候 变 化 
HRS (Copernicus Climate Change Service, C3S) 提供 
的 全 球 300 m 分 辨 率 土 地 利用 产品 (CCI-LC 与 
C3S-LC)。 按 照 李 元 春 等 "的 土地 覆 被 分 类 方法 ， 
将 土地 覆 被 重 分 类 为 草地 、 林 地、 裸 地 水 体 、 永 久 
冰雪 、 人 类 占用 地 共 6 类 ,其 中 人 类 占用 地 包括 耕地 
与 城区 。 
1.3 研究 方法 

根据 夏 自强 等 “在 该 流域 的 研究 ,将 四 季 划 分 


为 : 春 (4 一 5 月 )、 夏 (6 一 8 月 )、 秋 (9 一 10 月 )、 冬 (11 
月 一 次 年 3 月 )。 首 先 ,利用 伊 巴 流域 不 同时 间 尺 度 
的 ET 数据 ,计算 面 平 均值 得 到 上 中、 下 游 及 全 流域 
的 ET 值 ; 然 后 借助 Mann-Kendall(M-K) 法 于 ,分 别 
计算 逐 栅 格 上 .中 .下游 及 全 流域 的 蒸 散发 的 Z 值 ; 
基于 Theil-Sen median 趋势 分 析 呈 计算 上 中、 下 游 
月 尺度 ET 与 年 尺度 环境 要 素 的 年 际 变 化 率 。 最 后 
计算 每 个 栅 格 2010—2020 年 相对 2000—2010 年 多 
年 平均 年 ET 的 变化 率 ,分 析 不 同年 代 间 ET 的 空间 
变化 。 

根据 流域 2000 年 和 2020 年 的 土地 利用 数据 ， 
依据 土地 利用 转移 矩阵 原理 "制作 流域 土地 利用 
变化 图 ,分 析 土 地 覆 被 变化 对 流域 ET 的 影响 。 为 了 
确定 控制 不 同 区 域 ET 的 主要 环境 因素 ,基于 Pear- 
son 相 关系 数 法 "进行 月 ET 与 NDVI. 气 温 、 风 速 . 相 
对 湿度 、 土 壤 湿 度 5 个 环境 影响 因素 之 间 的 相关 性 
分 析 。Pearson 相关 系数 + 可 以 衡量 两 变量 间 线 性 相 
关 性 ,Irle[0,11 ,lr 二 0.8 时 为 高 度 相关 ,0.5<Irl<0.8 为 
中 度 相 关 ,0.3<lrl<0.5 为 低 度 相 关 ,Irl<0.3 为 不 相 
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关 。 根 据 年 尺度 水 量 平衡 ”” ,研究 伊 巴 流域 生态 
系统 水 资源 供给 量 的 变化 。 


2 结果 与 分 析 


21 流域 蒸 散 发 时 空 变 化 规律 分 析 
2.1.1 年 尺度 

(1) 多 年 平均 ET 的 空间 分 布 

流域 多 年 平均 ET 具有 明显 的 空间 分 布 特征 (图 
3) ,东南 高 .西北 低 ,最 大 值 可 达 1028.6 mm, 而 最 小 
值 仅 为 117.6 mm。ET 高 值 区 主要 位 于 卡拉 塔 尔 河 、 
阿克苏 河上 游 .阿拉 木 图 南部 、 喀 什 河 和 巩 乃 斯 河 
流域 等 植被 覆盖 度 较 高 的 区 域 ,其 土地 利用 类 型 多 
为 林地 或 耕地 ;ET 低 值 区 主要 分 布 在 巴 湖 周边 平原 
地 区 及 伊犁 河中 游 低 海拔 地 区 。 
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图 3 2000 一 2020 多 年 平均 年 ET 的 空间 分 布 特征 
Fig. 3 Spatial distribution characteristics of multi-year 


average evapotranspiration from 2000 to 2020 


全 流域 多 年 平均 年 ET 为 249.7 mm, MEP, 
下 游 多 年 平均 年 ET 则 分 别 为 439.0 mm 317.9 mm, 
201.1 mm, 比例 近 似 于 4:3:2, 上 中、 下 游 多 年 平均 
年 ET 存在 显著 差异 。 

(2) ET 变化 趋势 空间 分 布 

从 图 4 来 看 , 伊 巴 流域 大 部 分 区 域 年 ET 呈 增 加 
趋势 ,其 中 在 上 游 伊犁 河谷 及 天 山 山脉 的 大 部 分 地 
区 通过 99% 的 显著 性 检验 ;在 伊犁 河谷 的 周边 区 
域 ,年 莹 散发 量 呈 下 降 趋势 ,但 不 显著 ;上 游 蒸 散发 
的 变化 趋势 与 王 倩 等 ”提出 的 上 游 植被 生长 季 ND- 
VI 的 空间 变化 特征 具有 一 致 性 。 在 伊犁 河中 游 地 
区 , 绝 大 部 分 面积 的 年 ET 变化 趋势 不 显著 ;在 伊犁 
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图 4 2000 一 2020 年 蒸 散 发 的 M-K 趋 势 检验 Z 值 的 空间 分 布 
Fig.4 Spatial distribution of M-K trend test Z-values for 
evapotranspiration from 2000 to 2020 


河 下 游 的 三 角 洲 周边 地 区 及 巴 湖 的 东北 部 ,年 ET 呈 
增加 趋 热 , 且 通 过 了 95% 的 显著 性 检验 。 

(3) 上 上、 中、 下 游 年 ET 的 变化 趋势 

伊 巴 流域 上 中、 下 游 和 全 流域 年 ET 的 年 际 变 
化 (图 5) 规 律 相 似 , 均 呈 增 加 趋势 ,其 在 2000 一 2020 
年 间 的 平均 增幅 分 别 为 18.32 mm (10a)”',8.07 mm’ 
(10a) ,10.35 mm: (10a) 4112.22 mm*(10a)"。 上 、 
下 游 和 全 流域 年 ET 的 M-K 趋势 检验 Z 值 分 别 为 
1.74、1.96 和 1.96, 分 别 通过 了 90%、95% 和 95% 的 显 
著 性 检验 ,中游 蒸 散 发 量 增加 趋势 不 显著 。 
2.1.2 季节 尺度 

(1) 季节 尺度 ET 的 空间 分 布 

由 图 6 可 知 , 全 流域 春 . 夏 秋冬 四 季 的 多 年 平 


均 日 蒸 散发 量 分 别 为 :0.91 mm、0.92 mm .0.60 mm 
和 0.72 mm, ,而 上 游 四 季 的 分 别 为 :1.36 mm 2.47 


540 - 一 上游 —— 中 游 —— Pe 全 流域 


蒸 散 发 量 /mm 


2000 2002 2004 2006 2008 2010 2012 2014 2016 2018 2020 


年 份 


图 5 流域 上 .中 .下游 年 薰 散发 总 量 年 际 变化 
Fig. 5 Interannual variation of total evapotranspiration in the 


upper, middle and lower basin 
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图 6 2000 一 2020 年 四 季 蒜 散发 空间 分 布 特征 


Fig.6 Spatial distribution characteristics of seasonal evapotranspiration from 2000 to 2020 


mm 0.98 mm 和 0.89 mm, 中 游 的 分 别 为 :1.12 mm, 
1.37 mm .0.68 mm 和 0.85 mm, 下 游 的 分 别 为 :0.79 
mm 0.54 mm,0.51 mm 和 0.67 mm。 就 全 流域 而 言 ， 
春 、 夏 两 季 的 多 年 平均 日 蒸 散发 量 儿 乎 相同 ,秋冬 
两 季 的 则 小 一 些 。 上 ,中游 均 是 夏季 的 多 年 平均 日 
蒸 散 发 量 最 大 ,而 下 游 则 是 春季 的 多 年 平均 日 蒸 散 
发 量 最 大 。 从 高 程 和 NDVI 的 空间 分 布 (图 2a, 图 
2c) 来 看 ,四 季 ET 高 值 区 均 主 要 分 布 在 海拔 较 高 且 
植被 覆盖 度 较 高 的 地 区 ,但 在 夏 、 秋 两 季 ,伊犁 河谷 
地 区 也 呈现 较 高 的 日 均 蒸 散发 量 , 这 是 由 于 伊犁 河 
谷地 区 的 土地 利用 类 型 主要 为 耕地 , 且 当 地 在 3 月 
下 旬 至 4 月 上 旬 才 进入 播种 阶段 关 5 , 故 春 季 伊 犁 
河谷 的 日 均 蒸 散发 量 较 低 ,而 伊犁 河谷 外 天 然 植被 
履 盖 度 较 高 的 区 域 , 则 在 春季 呈现 出 较 高 的 日 均 蒸 
散发 量 。 

(2) 季节 尺度 蒸 散 发 变化 趋势 空间 分 布 

2000—2020 年 四 季 阁 散发 M-K 趋 势 检 验 Z 值 
的 空间 分 布 见 图 7。 上 游 伊犁 河谷 地 区 蔡 散 发 在 
夏秋 两 季 呈 增加 趋势 ,天 山 山脉 地 区 落 散 发 在 春 、 


夏秋 三 季 旦 增加 趋势 , 且 两 区 域 的 大 部 分 面积 通 
过 了 99% 的 显著 性 检验 ;中 游 大 部 分 地 区 蒸 散发 在 
四 季 中 均 呈 增加 趋势 ,但 不 显著 ;下 游 巴 湖 周边 区 
域 的 蒸 散发 在 春 E .秋季 呈 增 加 趋势 ,日 大 部 分 面 
积 通过 了 95% 的 显著 性 检验 。 冬 季 流 域 上 中 游 .中 
下 游 交 界 处 的 大 片区 域 蒸 散发 旺 不 显著 上 升 趋势 ， 
而 巴 湖北 部 局 部 地 区 的 蒸 散发 旺 下 降 趋 势 , 日 通过 
了 95% 的 显著 性 检验 。 
2.1.3 月 尺度 由 图 8 可 知 , 上 、 中 下游 ET 的 年 内 
分 配 呈 现 不 同 变化 特征 ,上 中游 ET 的 年 内 分 配 为 
单 峰 ,其 中 上 游 在 7 月 达到 峰值 81.1 mm, 中游 在 6 
月 达到 峰值 43.5 mm。 上 游 ET 的 峰值 相 较 于 中 游 
有 偏 高 偏 晚 的 分 布 特征 ,与 上 中游 作物 种 植 结构 
和 种 类 不 同 有 关 "”""。 下 游 ET 的 年 内 分 配 为 双 峰 ,在 
3 月 达到 峰值 后 开始 下 降 至 9 月 ,而 后 又 开始 上 升 至 
11 月 ,达到 男 1 个 峰值 后 再 次 下 降 , 其 中 3 月 和 11 月 
的 峰值 分 别 为 29.9 mm 和 24.4 mmo 

流域 不 同月 份 蒸 散发 的 年 际 变化 率 见 表 1。 流 
域 上 、 下 游 的 年 ET 旦 上 升 趋势 , 主 要 是 由 于 春 、 夏 、 
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Fig. 7 Spatial distribution of the M-K trend test Z-values for the seasonal evapotranspiration from 2000 to 2020 


2.1.4 ET 年 代 际 空间 变化 规律 ”由 图 9 可 知 ,在 上 


90[ A 一 上游 一 一 中 游 一 下 游 R 、 y 、 、 
al 游 伊犁 河谷 地 区 和 下 游 三 角 洲 地 区 ET 出 现 了 显著 增 
EKG 加 ,增幅 超过 了 10%。 三 角 洲 周边 区 域 ET 增幅 超过 
Z so 5%, 这 一 变化 趋势 与 图 4 的 M-k 的 空间 分 布 一 致 。 在 
Hi EERE KA AOK PERE IBA Be EIRATA 
20| 河谷 的 局 部 区 域 ,年 平均 ET 则 出 现 了 明显 的 下 降 , 降 
下 幅 超过 5%。 其 中 卡 普 丛 盖 水 库 南 部 ET 下 降幅 度 超过 


了 5%, 其 原因 可 能 是 2005 年 后 水 库 周边 由 ET 值 较 大 
的 耕地 持续 向 ET 值 较 小 的 天 然 植被 转化 ”” 。 
2.2 流域 蒸 散 发 影响 因子 变化 特征 

由 图 10 可 知 , 巴 湖 的 北部 和 伊犁 三 角 洲 周边 有 
大 片 的 裸 地 变 为 草地 ,而 在 三 角 洲 的 内 部 则 存在 森 
林 退 化 为 草地 。 在 上 游 较 为 明显 的 土地 利用 转换 


月 份 
图 8 2000 一 2020 年 多 年 平均 蒸 散发 的 年 内 分 配 


Fig. 8 Intra-annual distribution of multi-year average 


evapotranspiration from 2000 to 2020 


#k(4—10 H ) 季 ET 的 年 际 变 化 率 大 于 0, 其 中 上 游 9 


月 及 下 游 4 月 和 9 月 的 ET 变化 率 分 别 通过 95% 和 
99% 的 显著 性 检验 。 与 其 他 3 个 季节 相 比 ,冬季 (11 
月 一 次 年 3 月 ) 多 数 月 份 ET 的 年 际 变 化 率 小 于 0 且 
值 较 小 , 故 对 年 ET 的 变化 趋势 影响 不 大 。 


为 草地 变 为 林地 ,主要 分 布 在 天 山 山脉 处 ,一 定 程 
度 上 解释 了 该 地 区 燕 散发 的 上 升 ,此 外 ,草地 变 为 
人 类 占用 地 也 是 较为 明显 ,这 一 点 与 刘 婉 如 ”和 王 
宏伟 等 ”的 研究 成 果 一 致 。 
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表 1 2000—2020 年 月 蒸 散 发 的 年 际 变化 率 


Tab.1 Interannual trends in monthly evapotranspiration from 2000 to 2020 /mm'a) 
1 月 2 月 3 月 4 月 5 月 6 月 7 月 8 月 9 月 10 月 11 月 12 月 
上 游 0.011 -0.041 -0.038 0.264 0.191 0.239 0.455 0.249 0.25° 0.151 -0.039 0.020 
中 游 0.047 -0.036 0.040 0.247 0.272 -0.023 0.008 0.094 0.224 0.052 -0.033 0.019 
下 游 -0.112 -0.150 0.055 0.425" 0.178 0.284 0.245 0.148 0.267" 0.196 -0.086 -0.022 
注 :*# 通 过 99% 显 著 性 检验 ;* 通 过 95% 显 车 性 检验 。 下 同 。 
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图 9 2010 一 2020 年 相对 2000 一 2010 年 注 : 图 中 仅 显示 面积 变化 >500 km 的 土地 利用 变化 类 型 。 


多 年 平均 年 ET 变化 率 
Fig.9 2010-2020 relative to 2000-2010 multi-year average 


annual ET rate of change 


由 表 2 可 知 , 全 流域 的 风速 在 21 a 内 变化 趋势 
不 明显 ,而 降水 量 与 气温 在 21 a 里 呈 下 降 趋势 ,其 
中 ,中 游 降水 降幅 为 5.77 mm-a, Hil it 95% AY ib 
著 性 检验 。 全 流域 的 相对 湿度 .土壤 湿度 和 NDVI 
均 呈 增加 趋势 ,其 中 下 游 相 对 湿度 与 上 下游 的 ND- 
VI 增幅 均 通 过 95% 的 显著 性 检验 ,日 NDVI 的 变化 
趋势 与 段 伟 利 等 研究 中 得 出 的 2000 年 后 流域 内 
荒地 面积 减少 的 趋势 一 致 。 降 水 和 土壤 湿度 相反 
的 变化 趋势 反映 出 可 能 是 人 类 的 引水 灌溉 活动 造 
成 土壤 湿度 的 增 大 ,这 有 助 于 干旱 区 植被 的 生长 ， 
且 在 一 定 程度 上 也 促进 了 流域 ET 的 增 大 。 从 图 8、 
图 11 可 以 看 出 ,上 中游 ET 的 年 内 分 配 与 气温 和 
NDVI 的 年 内 分 配 相似 ,下 游 ET 的 年 内 分 配 则 与 土 


图 10 土地 利用 转移 空间 分 布 


Fig. 10 Spatial distribution of land use transfer 


壤 湿 度 和 相对 湿度 的 年 内 分 配 相 似 , 且 上 .中 游 多 
山区 ,植被 覆盖 度 较 高 ,下 游 多 平原 旦 植被 覆盖 度 
较 低 , 杀 异 的 气候 条 件 和 地 形 地 貌 可 能 使 上 、 中 、 下 
游 的 ET 受 不 同 因素 影响 。 
2.3 不 同 区 域 ET 的 影响 因素 

MOD16 的 ET 产品 是 基于 Penman-Monteith(PM) 
公式 "计算 ,运用 PM 公式 需要 气象 .植被 .土壤 水 
分 等 条 件 ”。 学 者 们 研究 了 基于 PM 公式 的 实际 
蒸 散 发 的 影响 因素 ,发 现 不 同 地 区 ET 的 主要 影响 因 
素 不 同 汪 2。 为 明确 影响 本 流域 内 ET 变化 的 主要 
因素 ,本文 基于 Pearson 相关 系数 法 ,计算 月 ET 与 5 
个 环境 影响 要 素 的 相关 系数 ,得 到 相关 系数 的 空间 
分 布 ( 表 3、 图 12)。 在 流域 的 上 中游,NDVI 和 气温 
与 ET 达到 中 等 程度 正 相 关 , 并 通过 了 99% 的 显著 性 


表 2 2000—2020 年 间 环 境 因子 的 年 际 变化 率 


Tab.2 Interannual trends of environmental factors from 2000 to 2020 


区 域 降水 量 /(mma) 相对 湿度 /(%…a') S/C +a!) 风速 /ms …a') 土壤 湿度 /kg'm > at) ~~ NDVI/(10a)" 
上 游 -1.65 0.06 -0.013 -0.002 0.27 0.016' 
中 游 -5.77 0.09 -0.016 -0.003 0.10 0.010 
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图 11 2000 一 2020 年 多 年 平均 环境 影响 因素 的 年 内 分 配 


Fig. 11 Intra-year distribution of multi-year average environmental impact factors from 2000 to 2020 


检验 ,表明 上 .中 游 的 ET 可 能 主要 受 这 2 个 因素 影 
响 , 且 会 随 着 这 2 个 因素 的 增 大 而 增 大 。 具 体 到 特 
定 区 域 来 看 ,伊犁 河谷 地 区 蒸 散发 与 气温 呈 显 著 的 
正 相 关 , 而 天 山 山脉 地 区 ET 则 同时 与 NDVI 和 气温 
均 呈 正 相 关 ( 图 12a, 图 12b)。 在 下 游 ET 的 影响 因 
素 中 , 仅 有 土壤 湿度 与 ET 达到 低 度 正 相 关 , 故 推测 
下 游 ET 主 要 受 土 壤 湿 度 控 制 , 即 一 方面 土壤 水 分 增 
加 改善 植被 生长 条 件 促 进 植被 蒸腾 , 另 一 方面 也 会 
促进 土壤 自身 蒸发 ;根据 PM 公 式 , 在 其 他 条 件 不 变 
的 情况 下 , 当 土 壤 水 分 在 凋 萎 系数 和 饱和 含水 量 之 
间 增 加 时 , 蔡 散 发 所 受 的 胁迫 作用 将 减轻 ,土壤 水 
分 胁迫 系数 增 大 进而 促使 ET 增 大 中 。 具 体 到 特定 
区 域 来 看 ,下 游 巴 湖北 岸 和 南岸 的 土壤 湿度 与 ET 的 
相关 性 分 别 达 到 了 低 度 和 中 度 相关 ,而 在 伊犁 河 三 


表 3 流域 上 ,中 .下游 蒸 散发 量 与 环境 因子 的 相关 性 
Tab.3 Correlation between evapotranspiration 
and environmental factors in the upper, middle 


and lower basin 


影响 因素 上 游 中 游 下 游 
NDVI 0.735" 0.681" 0.012 
风速 -0.283" -0.033 0.129” 
温度 0.790” 0.552” -0.129° 
相对 湿度 -0.288” -0.095 0.160° 
土壤 湿度 -0.246" -0.094 0.304 


角 洲 的 湿地 沼泽 区 和 阿克苏 河 与 卡拉 塔 尔 河 上 游 
的 林地 耕地 区 域 (图 2b) ,土壤 湿度 与 ET 呈 负 相关 
或 不 相关 关系 ,类 似 的 现象 也 体现 在 相对 湿度 与 ET 
的 相关 系数 空间 分 布 上 。 出 现 这 种 现象 的 原因 可 
能 是 由 于 该 区 域 较 大 的 降雨 量 * 或 人 为 灌溉 活动 
使 得 芋 散 发 所 需 的 供水 条 件 得 到 保障 ,供水 条 件 
(土壤 湿度 ) 不 再 是 决定 这 些 地 区 落 散 发 的 关键 要 
素 ,而 植被 覆盖 度 和 气温 才 是 影响 ET 大 小 的 主要 因 
素 。 风 速 与 ET 的 相关 系数 在 大 部 分 区 域 未 通过 
95% 的 显著 性 检验 , 仪 在 局 部 地 区 与 菩 散 发 呈现 显 
著 的 相关 关系 (图 12c) ,说 明 在 多 数 地 区 ,风速 并 不 
是 ET 的 主要 影响 因素 。 


3 讨论 


梁 红 闪 等 研究 了 仇 巴 流域 2000 一 2014 年 ET 
的 变化 特征 ,认为 流域 总 体 呈 下 降 趋 势 , 旦 未 来 以 
持续 性 减少 为 主 。 而 根据 本 文 研究 结果 ,全 流域 年 
ET 的 年 际 变化 呈 增 加 趋势 ,2000 一 2020 年 间 的 平均 
增幅 为 12.22 mm.(10aj。 图 5 说 明了 2014 年 后 , 流 
域 上 中、 下 游 年 蒸 散发 量 上 升 趋势 明显 ,经 计算 
E 、 中、 下 游 2015—2020 年 年 均值 分 别 为 466.59 
mm, 334.26 mm, 217.38 mm, 均 显 著 大 于 2000 一 
2014 年 年 平均 值 (427.91 mm 311.36 mm, 194.58 
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注 : 仅 显示 通过 95% 显 著 性 检验 的 数据 。 
图 12 蒸 散发 与 环境 因子 的 相关 系数 的 空间 分 布 


Fig. 12 Spatial distribution of correlation coefficients between evapotranspiration and environmental factors 


mm)。FET 的 持续 升 高 会 直接 影响 到 伊 巴 流域 生态 
系统 的 水 资源 供给 服务 能 力 , 故 基 于 水 量 平 衡 原 
理 , 可 以 计算 出 流域 生态 系统 水 资源 供给 服务 量 在 
ET 显著 增加 背景 下 的 变化 趋势 (图 13). 2000 年 
来 ,流域 的 年 降水 量 呈 下 降 趋 势 ( 表 2) ,而 流域 的 年 
ET 呈 增 加 趋势 (图 4) ,这 导致 流域 的 生态 系统 水 资 
源 供给 服务 总 量 出 现下 降 的 趋势 ,其 中 中 游 的 下 降 
趋势 最 明显 。 此 外 ,流域 上 游 的 生态 系统 水 资源 供 
给 服务 总 量 在 一 些 年 份 已 经 出 现 了 负 值 ,到 2020 年 
这 一 缺口 甚至 达到 了 100 mm ,而 中 .下 游 的 生态 系 
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图 13 2000 一 2020 年 生态 系统 的 水 资源 
供给 服务 量 的 年 际 变化 
Fig. 13 JInterannual variation in water supply services of 
ecosystems from 2000 to 2020 


统 水 资源 供给 服务 总 量 也 在 2020 年 出 现 了 负 值 。 
近年 来 上 游 伊犁 河谷 地 区 的 土地 利用 类 型 部 分 章 
地 向 耕地 转变 ” ,在 一 定 程度 上 加 剧 了 水 资源 的 
消耗 。 中 游 地 区 随 着 卡 普 恰 盖 水 库 的 建立 ,周边 区 
域 的 草地 等 几乎 全 部 逐渐 转变 为 农业 用 地 ,加 剧 了 
中 游 地 区 水 资源 的 消耗 ,同时 对 维持 巴 湖水 位 也 会 
造成 不 利 影 响 ”;1991 年 后 ,中 、 下 游 的 许多 耕地 面 
积 处 于 押 荡 状态 ,后 经 过 一 系列 土地 改革 措施 *”， 
耕地 面积 开始 大 大 增加 ,与 此 同时 ,下 游 巴 湖 附近 
的 许多 裸 地 也 开始 向 草地 转变 ,土地 覆 被 的 变化 也 
BRENT ARATE SEK 


4 结论 


为 了 研究 伊 巴 流域 ET 的 演变 规律 及 其 环境 影 
响 因素 ,利用 MOD16 产 品 结合 NDVI 气象 数据 和 土 
壤 湿 度 等 数据 得 到 以 下 结论 : 

(1) 伊 巴 流域 2000 一 2020 年 年 ET 在 空间 上 时 
现 自 上 游 向 下 游 逐 渐 减 少 的 分 布 ,上 中、 下游 多 年 
平均 年 ET 分 别 为 439.0 mm 317.9 mm 201.1 mm, Lt 
例 近 似 4:3:2。 其 中 上 PERFI H RAA E 
大 ,而 下 游 在 春季 最 大 。 流 域 上 中 游 蒸 散 发 量 的 年 
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内 分 配 均 为 “ 单 峰 型 ” ,峰值 分 别 在 7 月 与 6 月 ,下 游 
为 “ 双 峰 型 ” ,峰值 分 别 在 3 月 与 11 月 。 

(2) 下 游 的 年 ET 均 呈 显著 上 升 趋势 ,显著 区 域 
主要 分 布 在 上 游 的 伊犁 河谷 天山 山脉 与 下 游 的 伊 
Am = aae Ei 2010—2020 
年 
ww 

(3) 中 游 蒸 散发 与 气温 和 NDVI 呈 较 高 的 正 相 
关 性 ;流域 下 游 蒸 散 发 与 土壤 湿度 呈 低 度 正 相关 。 

(4) 受到 实际 蒸 散发 的 显著 增加 和 降雨 量 减少 
的 影响 ,流域 生态 系统 水 资源 供给 服务 总 量 在 
2000 一 2020 年 有 所 下 降 ,其 中 上 游 已 多 次 出 现 缺 
口 ,中 、 下 游 在 2020 年 开始 出 现 缺 口 ,急需 通过 控制 
用 水 总 量 .提升 用 水 效率 保障 水 资源 供需 平衡 。 
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Spatial and temporal variability of evapotranspiration and influencing factors 
in the Ili River-Balkhash Lake Basin 


YAO Jia, CHEN Qihui', LI Qiongfang'*, CUI Gang', ZHANG Liangjing’ 
(1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, Jiangsu, China; 
2. Yangtze Institute for Conservation and Green Development, Nanjing 210098, Jiangsu, China; 


3. College of Environment, Hohai University, Nanjing 210098, Jiangsu) 


Abstract: From 2000 to 2020, remote sensing data, such as MODIS and FLDAS, the Mann- Kendall method, 
Theil-Sen median trend analysis, Pearson correlation analysis, and the water balance principle, were applied to 
explore the spatial distribution and temporal variation patterns of actual evapotranspiration in the Ili River- 
Balkhash Lake basin, along with its main influencing factors. In addition, changes in the ecosystem water supply 
were discussed. The results were as follows: (1) the upper, middle, and lower Ili River-Balkhash Lake basins had 
an average annual evapotranspiration of 439.0 mm, 317.9 mm, and 201.1 mm, respectively. The daily 
evapotranspiration in the upper and middle basins was the greatest in the summer, while the lower basin had the 
greatest daily evapotranspiration in the spring. The intra- annual distribution of evapotranspiration in the upper 
and middle basins was “unimodal,” with peaks in July and June, respectively. It was “bimodal” in the lower 
basin, with peaks in March and November, respectively. (2) The annual evapotranspiration increased in both the 
upper and lower basins, with the highest concentrations in the upper Ili River valley, Tianshan Mountains, and 
near the Ili River delta in the lower basin. Compared to the 2000s, the average annual evapotranspiration in the Ili 
River valley and Ili River delta increased by more than 10% in the 2010s. (3) Based on the Pearson correlation 
analysis, the upper and middle basins showed high positive correlations with temperature and NDVI, whereas the 
lower basin showed high positive correlations with soil moisture. (4) From 2000 to 2020, the basin's total water 
supply service decreased, with several shortfalls in the upper basin and shortfalls in the middle and lower basins 
beginning to appear in 2020. The balance of water supply and demand must be ensured by controlling the total 
amount of water being used and by improving water use efficiency. 

Keywords: Ili River-Balkhash Lake Basin; evapotranspiration; temporal and spatial changes; driving factors; 


Pearson correlation analysis. 


